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ABSTRACT 

We report on the detection with the HIFI instrument on board the Herschel satellite of the two hydrogen chloride isotopologues, H^'Cl and 
tf'Cl, towards the massive star-forming region W3 A. The 7=1-0 line of both species was observed with receiver lb of the HIFI instrument at 
~625.9 and ~624.9 GHz. The different hyperfine components were resolved. The observations were modeled with a non-local, non-LTE radiative 
transfer model that includes hyperfine line overlap and radiative pumping by dust. Both effects are found to play an important role in the emerging 
intensity from the different hyperfine components. The inferred H'^'Cl column density (a few times ~10'* cm"-), and fractional abundance relative 
to H nuclei (~7.5xlO""'), supports an upper limit to the gas phase chlorine depletion of r;200. Our best-fit model estimate of the H^^Cl/H^'Cl 
abundance ratio is a;2. 1±0.5, slightly lower, but still compatible with the solar isotopic abundance ratio (~3.1). Since both species were observed 
simultaneously, this is the first accurate estimation of the ['^C1]/[^^^C1] isotopic ratio in molecular clouds. Our models indicate that even for large 
line opacities and possible hyperfine intensity anomalies, the H'^'Cl and H'^'Cl 7=1-0 integrated line-intensity ratio provides a good estimate of the 
'^Cl/^^^Cl isotopic abundance ratio. 
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1. Introduction 

. Chlorine has two stable isotopes (^^Cl and ^^Cl) and an ioniza- 
' tion potential of 12.97 eV (i.e., slightly below that of hydrogen). 
Hence, it can be ionized by UV photons (912-956 A) in diffuse 
clouds and in the edges of photon-dissociation regions (PDRs). 
Once ionized, Cl^ reacts with molecular hydrogen exothermi- 
. cally to form HCl^, a process that initiates the chemical reac- 
■ tions of chlorine. In cloud interiors, HCl^ can be formed by re- 
actions starting with neutral CI and H3 . The chemistry of chlo- 
rine in interstellar clouds has been the subject of various works 
(see Neufeld & Wolfire 2009). These studies predict that hydro- 
gen chloride (HCl) is the most abundant Cl-bearing molecule in 
dense clouds. 
] The HCl hyperfine lines (see Sect. 3.2) can be resolved 
. in interstellar sources only with heterodyne receivers equipped 
' with high spectral resolution spectrometers. The HCl 7=1-0 
' line at ~625.9GHz was first detected with the Kuiper Airbone 
' Observatory towards Orion (Blake et al., 1985) and followed 
' by detections in S gr B2 dZ muidzina set al.L Il995h . several po- 
sitions in OMC-1 dSchifkeet al.. .199lh . and in Mon R2 using 
the Caltech Submillimeter Observatory with good atmospheric 
transparency (Salez et al. 1996; SFL96 hereafter). The inferred 
HCl column densities are in the range lO'^-lO'^^cm"^. SFL96 
also presented the first detection of H^'^Cl towards Orion. Since 
the chemical reactions involving HCl are relatively well under- 
stood (see Sect. [T2l l. the H^^Cl/H^^Cl abundance ratio should 
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provide a good measure of the ■'^Cl/^^Cl isotopic ratio. Both ^^Cl 
and ^^Cl nuclei are believed to form in the last burning stages of 
massive stars (>10 M©) and by means of "explosive nucleosyn- 
thesis" during supernovae detonation (e.g., Woosley & Weaver 
1995). Therefore, observations of H^^Cl and H^'^Cl, and accurate 
measurements of the •'^Cl/^^Cl ratio in different environments, 
can provide some insight into the chemical evolution of both 
isotopes, thus into Galactic chemical evolution. The species HCl 
was also detect ed recently towards the carbon-rich evolved star 
IRC+10216 bv lCernicharoeTaLldMol) . 

Using HIFI, the He terodyne Instrument for the Far- 
Infrared (Ide Graauw et all 12010). on board the Herschel Space 
Observatory (Pilb ratt et al.L 12010), we present in this Letter the 
detection of the 7=1-0 rotational transition of H^^Cl and 
H^^Cl towards the massive star- forming region W3. The broad 
frequency coverage of HIFI allows us to observe both isotopo- 
logues of HCl with the same relative calibration in a wide variety 
of astronomical environments. Here we present the first accurate 
determination of the ^^C\f^^C\ isotopic ratio by a detailed model 
of the excitation of the hyperfine levels developed using a non- 
local radiative transfer code. The effect of line overlaps between 
the hyperfine components, and radiative pumping by dust pho- 
tons are discussed and modeled in detail. 



2. Observations and data reduction 

All spectra presented here were taken during the performance 
verification (PV) phase of HIFI (de Graauw et al. 2010). Both 
H^'^'Cl and H^'^Cl 7=1-0 lines were observed in the Band lb re- 
ceiver using the Wide Band Spectrometer (WBS), which pro- 
vides ~4GHz of bandwidth and ~1.1 MHz of channel spectral 
resolution (or a velocity resolution of ~0.5 km s"' at ~626 GHz). 



J. Cemicharo, J.R. Goicoechea, F. Daniel et al.: The ""'Cl/'^Cl ratio in molecular clouds 



'I'a 


_ ' 1 


0.2 


J=l-0 


0.1 


7 


" il 





WW 


ft 


k 




3f li\ 


\r 


' 1 







J=l-0- 









625 




y(GHz) 


Ta 


: 


^7 

H Cl 


0.2 


- 


- 


0.1 


- 


Pf 


- 




[ I 


J 


If 


^4, [ 





r 


^ 




1 V/\}, 




- 


1 1 1 1 1 1 1 1 1 1 1 1 H 



625.5 








626 


0.2 


■ i 


1 


35' ' ' - 
H Cl - 


0.1 


- 


u 




- 





H 


r 




1 


IVI 




- 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 - 



624.9 624.95 625 
i/(GHz) 



625.05 



625.85 625.9 625.95 
y(GHz) 



626 



Fig. 1. Detection of tf ^Cl and H"C1 7=1-0 lines towards W3 A 
Hii region. Arrows shows the relative line strength of each HFS 
component. The length of the arrows are proportional to the ex- 
pected intensities in the LTE optically thin limit. Both lines were 
observed simultaneously in a line survey of band lb of HIFI, 
hence, have the same calibration accuracy. Spectral resolution 
was smoothed to 2 MHz (^1 kms"'). 
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Fig.2. H^^Cl, H"C1 7=1-0, and C'^0 7=5-4 lines observed 
with HIFI towards W3 A. Continuous curves show closest fit- 
ting radiative-transfer model line profiles (see text). The spectral 
resolution is ^1 kms"' 



The telescope, which has a half -power beam-width (HPBW) of 
-35" at -626 GHz, was centered on azooo = 02''25'"43.51', 
^2000 = 62°06'13" (a position close to W3 A IRS2 and 2a). A 
complete spectral scan of the Band lb was taken at this position 
during PV phase and several lines from different molecules have 
been identified so far. The horizontal (H) and vertical (V) polar- 
ization receivers were averaged after rescaling the V one using 
the HCO+ J=7-6 line intensity at 624.21 GHz, i.e., close to that 
of HCl isotopologues. The data were first processed with HIRE 
software (Ott et al. 2010), and then exported to CLASS where 
standard data reduction routines were carried out. We checked 
that the target lines are not contaminated by lines from the other 
side band using the different frequency settings. The rms noise at 
~626GHz is ~30mK (antenna temperature) per 0.5 km s' res- 
olution channel. Hence, the H^^Cl and H^^Cl 7=1-0 lines are de- 
tected at lOcr and 6cr levels, respectively and with the same rel- 
ative calibration. At this level of sensitivity, only the most abun- 
dant species are detected. In particular, we detected 20 lines in 
the entire band lb above 3cr. The brightest is the CO 7=5-4 line, 
followed by the ground-state line of ortho-HaO and a few lines 
from formaldehyde, methanol, HCO^, and HCN. Hence, we are 
confident that the observed HCl line profiles are not blended with 
other spectra l features. By exam i ning our own and public spec- 
ti-al catalogs (Mu Uer et al.1 12001". 2005: Pickett et al.[ ll998h . we 
also checked that lines in the signal band from other molecules 
do not blend with the hyperfine components of both HCl iso- 
topologues. 

The total integration time was 2 minutes. Figures [T] and |2] 
show the resulting line profiles (data smoothed to a spectral res- 
olution of ^Ikms '). To compare with our models, the fol- 
lowing expression for the main beam efficiency was adopted, 
Imb - 0.72 exp(-(v/6)^) xO.96, where v is the frequency in THz, 
0.72 is rjinh for the telescope in the limit of frequency, and the 
factor 0.96 is the assumed forward efficiency of the telescope 
(M. Olberg priv. comm.). 



3. Results 

In terms of spectroscopy, the 7=3/2 nuclear spin of ^^Cl and of 
^^Cl splits the pure rotational transitions of H^^Cl and H-^^Cl into 
several hyperfine structure (HFS) components (see e.g., Cazzoli 
& Puzzarini 2004 and references therein). These hyperfine com- 
ponents are indicated as vertical arrows in Fig.[T] In the optically 
thin limit, these components follow a 2:3:1 intensity ratio (from 
the lowest to the highest frequency hyperfine component). It is 
clear from Fig. [1] that the observed ratios are close to 1:1:1 for 
the three hyperfine components of H^^'Cl. These ratios indicate 
that the H^^Cl hyperfine components are affected substantialy by 
opacity. Even for H^^Cl, the observed hyperfine line intensities 
( 1 : 1 :2) do not follow the expected ratios in an optically thin case. 
The intensity peak ratio of the strongest H^^'Cl and H^^Cl 
HFS components is ^1.5 and the integrated intensity ratio is 
^2+0.2. These val ues are lower than t h e sola r ^^Cl/^^Cl abun- 
dance ratio, ~3.1 (lAnders & Grevesso. Il989h . which suggests 
that optical depth effects could influence the observed HFS line 
intensity ratios (or that the H^^Cl/H-^^Cl abundance ratio is lower 
than the solar value). For completeness. Figure |2] compares the 
detected HCl lines with the C'^O 7=5-4 at 561.712GHz also 
observed with HIFI. The detection of C'^O, HCO+ and HCN 
mid-7 lines confirms the presence of warm and dense molecular 
gas towards the observed position. 



3.1. HCl excitation and ^^Cl/^^Cl abundance ratio 

The star-forming region W3 is located in the Perseus arm at a 
distance of 2.3 kpc and contains several young massive stars that 
ionize a natal molecular cloud creating Hn regions. In particu- 
lar, the near-IR sources IRS 2 and IRS 2a (OB stars) are believed 
to be the ionizing sources of the W3 A Hii region (Tieftrunk et 
al. 1995 and references therein). These sources ar e also the ori- 
gin of molecular outflows and are X-rays emitters dHofner et al.L 
I2002h . 
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To estimate the H^^Cl/H^'^Cl abundance ratio and analyze all 
possible effects affecting the emerging line profiles, we modeled 
the observed HCl and C'^O lines with our non-local and non- 
LTE r adiativ e transfer codes dO onzalez-Alfonso & Cernicharo, 
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119971; iGoicoechea et al.L ,2006; .Daniel & Cernichai'o, 
To take into account the blending of the HCl 7=1-0 hy- 
perfine components (i.e., to determine the opacity at each fre- 
quency when several lines overlap), we used the modelling ap- 
proach presented in Daniel & Cernicharo (2008) to interpret 
the HFS line emission from N2H+, HCN, and HNC. Rate co- 
efficients for the collisional excitation of HCl by He are taken 
from Neufeld & Green (1994), who also estimated the con- 
tribution to each specific hyperfine level. Like other light hy- 
drides that HIFI will observe, HCl has a large rotational constant 
(10.4 cm '), thus HCl rotational transitions have high sponta- 
neous radiative rates. The HCl critical densities are very high, 
ncr(J=l-0)-7xl0^cm"'', and only when n(H2)>ncr does colli- 
sional excitation dominate. This high value suggests that HCl 
line emission arises in dense molecular gas. Radiative pumping 
by dust photons may also be very important in determining the 
HCl level populations (see also SFL96) because of the increase 
in grain emissivity an d dust opac ity in the far-IR and submil- 
limeter domains ( Cernicharo et al.L 12006 a.b.) . Model predictions 
shown in Fig. |3] demonstrate that the inclusion of line overlap, 
dust pumping, and both effects together, modifies the relative in- 
tensity of each HFS component. We note in particular how the 
7=1-0 F -1/2-1/2 line (the HFS component with the weakest 
line strength and opacity) is enhanced with respect to the other 
components when radiative pumping is included. 

To reproduce the observed H^^Cl and H^^Cl line profiles 
and relative intensities, we assumed uniform physical conditions 
(nH:^10''cm-^ and T^^IOOK; taken from Helmich et al. 1994; 
Tieftrunk et al. 1995) and that the H^^Cl and H"C1 abundances 
are free parameters. The C'^O line was also analyzed with the 
same parameters to verify the model consistency. In particular, 
the C'^O 7=5-4 line was found to be optically thin (t^O.1), 
which allowed us to constrain the line-of-sight column density 
of material and also the line velocity dispersion (cr^l.3kms '). 
Gas and dust were assumed to coexist and be thermally coupled 
(Ti=T^). 

To reproduce the observed HCl line peak positions and their 
relative strengths, best fit solutions were obtained for an expand- 
ing shell of gas. The adopted velocity gradient (from 2.5 km s ' 
at the center to 0.5 km s"' at the edge) is consistent with the 
CO molecular outflows seen in the region (e.g., Hasegawa et al. 
1994). Optimal results were obtained for a H^^Cl column den- 
sity of a few times 10'"^ cm"^ (or an abundance of ~7.5xl0""' 
relative to total H). The inclusion of radiative pumping from 
dust (and line overlap to a lesser extent) allowed one to more 
accurately reproduce the observed HFS relative line intensities. 
Assuming that H^^Cl and C'^O arise from the same regions, we 
inferred a column density ratio of A^(H^5C1) - A^(C'^O)/20, by 
using the CO abundance determ ined for the region (~4xl0"^ per 
H nucleus; [Tielens et al.l (Il991h ) and assume a standard '^O/'^O 
isotopic ratio of 2600. Optimal fits were obtained by using an au- 
tomatic x^ procedure and are shown in Fig. |2] The H^^Cl/H^^Cl 
abundance ratio found in the models is ~2.1 (with a confidence 
interval within 1.6-3.1). Since H^^Cl lines are moderately opti- 
cally thick (t~12, 8, 4 for each HFS component), one does not 
expect the observed H^^Cl/H^'^Cl line intensity ratio to provide 
a direct measure of the H^^Cl/H^'^Cl abundance ratio. However, 
HCl critical densities are much higher than the gas density in 
most ISM clouds, and therefore the excitation temperature of the 
different HFS components remain proportional to the HCl col- 




Fig. 3. Results from non-local, non-LTE radiative transfer cal- 
culations for the H^^Cl 7=1-0 HFS components. The different 
curves show the effects of including line overlap and radiative 
pumping from dust photons in the emerging line intensities and 
relative HFS line ratios. 
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Fig. 4. Modeled HCl and H^^Cl 7=1-0 integrated fine intensity 
ratio (R) as a function of the assumed isotopic abundance ratio 
(see text). HCl abundances, indicated at the top left, are relative 
to H nuclei. Continuous lines correspond to models with a ve- 
locity gradient. Dashed lines correspond to a static model. 



umn density, even for optically thick lines. To conclude whether 
or not the line integrated intensity ratio is a good measure of the 
isotopic abundance ratio. Fig. 4 shows the modeled integrated 
intensity ratio as a function of the H^'^Cl/H^^Cl abundance ra- 
tio. Although the hyperfine components are optically thick in 
most models, the integrated line intensity ratio is proportional to 
the isotopic ratio for HCl abundances below 10"^. In Fig. 4, we 
also show the results for a static cloud. In this case, the opac- 
ities are larger but the integrated intensity ratio still provides a 
reliable measurement of the isotopic abundance ratio for HCl 
abundances below 10"'°. 



3.2. Chlorine chemistry in W3A Ha region 

In the UV-illuminated gas, the chlorine chemistry involves the 
reaction of Cl^ with molecular hydrogen to form HCL, which 
then reacts with H2 to produce H2CL. If the abundance of elec- 
trons is high, dissociative recombination of H2CL leads to the 
formation of HCl with a typically assumed branching ratio of 
~10%. In cloud interiors, atomic chlorine is mostly neutral, not 
ionized, so that the reaction of CI and H3 drives the formation 
of H2CL, which then reacts with CO and H2O leading to the 
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Fig. 5. Cloud-depth dependent photochemical models adapted 
to W3 A physical conditions. The gas density (n//=10^cm"'') 
and temperature (Ti=100K) are kept constant. The UV radia- 
tion field is x-^^"" ^nd the ionization rate due to cosmic rays 
is ^(H)=2.5xl0 '^s '. The predicted abundance of several Cl- 
bearing species are shown. A chlorine gas-phase abundance of 
1x10 "^ is used. The dashed line shows the expected HCl abun- 
dance with an undepleted abundance of [Cl] = 1.8xlO ''. 



formation of HCl. An alternative direct route to HCl in either 
hot gas or regions where vibrationally excited H2 is abundant 
{e.g., Agundez et al. 2010) is the CI -H H2 ^ HCl + H reaction 
which possesses an energy barrier of ~0.2 eV (iDobis & BensonL 
12002 ). The destruction of HCl is dominated both by pho- 
toionization and photodissociation and by reactions with C^ 
and H3 (the latter in the UV shielded gas). Previous observa- 
tional studies of HCl suggest that there has been a depletion 
of gas-phase chlorine in dense molecular clouds of a factor 
of >100 (SFL96) relative to the elemental chlorine abundance 
observ e d in diff'use clouds (JCl1^l .8xlO~^; ISavage & Sembach 
(Il996l) . ISonnentrucker et aP (l2006l) and references therein). To 
follow the HCl chemistry in the particular environment of W3 A 
Hii region and estimate the CI depletion we modeled the Cl- 
photochemistry using the Meudon PDR code (Le Petit et al. 
2006; Goicoechea & Le Bourlot 2007). The reaction network 
for Cl-bearing molecules includes the updated rates of Neufeld 
& Wolfire (2009). We note that the X-ray luminosity reported 
in the region (Lx~5xlO^' ergs"'; Hofner et al. 2002) is insuf- 
ficient for a "XDR-" rather than a "PDR-dominated" environ- 
ment. However, in cloud interiors X-ray photons may play an 
important role in the Cl-chemistry. 

Figure|5]shows the output of a model adapted to the physical 
conditions in W3 A. The UV radiation field produced by the OB 
stars in the region is simulated by an enhancement of 10'* times 
the mean interstellar radiation field (in Draine units). To repro- 
duce the inferred H^^Cl abundance (~7.5xl0 '*'), a gas-phase 
chlorine depletion of <200 is needed, HCl accounting for kTIWo 
of the CI nuclei in the gas phase. This conclusion is reached as- 
suming that the observed HCl arises in regions of large H2 col- 
umn density (Av>100; see Sect l3.1b . which is consistent with 
the submm continuum maps of the region (JafFe et al. 1983). If 
the observed HCl arises in regions of lower extinction, the CI 
depletion factor will obviously be lower. 



4. Conclusions 

We have presented the first detection of H-'^Cl and H^'Cl to- 
wards the W3 A Hii region. The inferred H^^Cl column density 
(a few times ~10''*cm"^) and fractional abundance (~7.5xl0"'" 
per H nucleus) provide an upper limit to the gas phase chlo- 
rine depletion of si200. This value is lower than that observed 
towards Orion hot core, but similar to that inferred towards 
Mon R2 (SFL96). Radiative transfer models including HFS line 
overlap and pumping by dust photons have been used to inter- 
pret the observations. The best-fit model provides a H^^Cl/H^^Cl 
abundance ratio of s;2.1, which is both lower than the so- 
lar value (~3.1) and lower than the previous estimate towards 
Orion (x!4-6; SFL96). On the other hand, it is similar to the 
[^^^C1]/["C1] ratio obtained in the IRC-H10216 circumstellar en- 
velope from [Na^^Cll/rNa^^Cn and fAl^^Cll/f Al^^Cll meas ure- 
ments dCernicharo & GueUnl [l987t ICernicharo et aLlboOOl) . 
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